Comparing and combining the Saskatoon, QMAP and COBE CMB maps 
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We present a method for comparing and combining maps with different resolutions and beam shapes, 
and apply it to the Saskatoon, QMAP and COBE/DMR data sets. Although the Saskatoon and 
QMAP maps detect signal at the 21<r and 40cr levels, respectively, their difference is consistent with 
pure noise, placing strong limits on possible systematic errors. In particular, we obtain quantitative 
upper limits on relative calibration and pointing errors. Splitting the combined data by frequency 
shows similar consistency between the Ka- and Q-bands, placing limits on foreground contamination. 
The visual agreement between the maps is equally striking. Our combined QMAP+Saskatoon map, 
nicknamed QMASK, is publicly available at www. hep. upenn.edu/^xuyz/qmask. html together with 
its 6495 x 6495 noise covariance matrix. This thoroughly tested data set covers a large enough 
area (648 square degrees — currently the largest degree-scale map available) to allow a statistical 
comparison with COBE/DMR, showing good agreement. 
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FIG. 1. The three maps that we will compare and combine are 
shown in Galactic coordinates. QMAP location in COBE map is 
shown in black. 



I. INTRODUCTION 

The cosmic microwave background (CMB) field is cur- 
rently enjoying a bonanza of new high-quality data jl]-[| , 
which has triggered a surge of new papers about the im- 
plications for cosmological parameters ]^-pT|. There is 
currently such wide interest in these cosmological results 



that it is tempting to temporarily ignore underlying as- 
sumptions. However, it is nonetheless important to con- 
solidate these gains by careful study of the many tech- 
nical analysis steps upon which they rest. This can be 
done at many levels. The experimental teams generally 
test for systematic errors at many steps in their data 
analysis pipeline, from data acquisition, cleaning, cali- 
bration and pointing reconstruction to mapmaking and 
power spectrum estimation. In addition, numerous de- 
tailed comparisons have been made between the angu- 
lar power spectra Ct measured by different experiments 
to determine whether they are all consistent 
However, such comparisons use only a very small frac- 
tion of the information at hand: average band powers, 
not the spatial phase information. A more powerful test 
(in the statistical sense of being more likely to discover 
systematic errors) involves a direct comparison of the sky 
maps from experiments that overlap in both spatial and 
angular coverage. 

Such a comparison is straightforward for maps with 
identical resolution and beam shape, simply testing 
whether the difference map is consistent with pure detec- 
tor noise. Such tests have been successfully performed for 
the COBE/DMR maps p7| . Unfortunately, comparisons 
are usually complicated by angular resolution differences 
between channels. Some experiments probe the sky in an 
even more complicated way, with, e.g., elliptical beams, 
double beams, triple beams, interferometric beams or 
complicated elongated software-modulated beams. Cor- 
related noise further complicates the problem. Despite 
these difficulties, accurate comparisons between different 
experiments are crucial. Some of the best testimony to 
the quality of CMB maps comes from the success of such 
comparisons in the past — between FIRS and DMR [18|, 
Tenerife and DMR Q, MS AM and Saskatoon pU20|, 
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two years of Python data |22| , three years of Saskatoon 
data [p3l , two flights of MS AM Hi and different channels 



of QMAP [E5-27|, Boomerang H and Maxima 



General methods have been developed for both com- 
paring ^Q,Eq] and combining |2^] arbitrarily complicated 
experiments. In this paper we will derive a technique for 
simplifying this task in practice, and apply it to com- 
pare and combine the three overlapping data sets shown 
in Figure § COBE/DMR, Saskatoon and QMAP. Our 
motivation is threefold: 

• To test and provide methods that can be used by 
experimental groups in the future. 

• To search for systematic problems that may be rel- 
evant to ongoing and future experiments. 

• To quality-test and publicly release the largest 
degree-scale map to date. 

We stress that combining maps is not just a matter of 
making pretty pictures. Power spectra from different ex- 
periments are routinely combined as if their sample vari- 
ance were independent. However, since this approxima- 
tion breaks down whenever the underlying maps overlap 
in spatial and angular coverage, the only correct way to 
compute their combined power spectrum is to extract it 
from the combined map. 

We present our methods in Section [fi|, present our re- 
sults in Section [II and summarize our conclusions in 



Section IV 



II. METHOD 

We use the methods for combining and comparing 
maps presented in |2^| , which are most easily expressed 
with matrix notation. Given two data sets represented 
by the vectors yi and y2, we write 



yi=Aix + ni, y 2 = A 2 x + n 2 . 



(1) 



Here the vector x contains the temperature of the true 
sky at various locations (pixels). Ai and A 2 are two 
known matrices incorporating the pointing strategy and 
beam shape of each experiment, ni and ri2 are two ran- 
dom noise vectors with zero mean and known covariance 
matrices Ni = (nin*) and N 2 = (n 2 ri2). It is convenient 
to define larger matrices and vectors 



A = 



Ai 



- / yi 



n 



A 2 J ' - \ y 2 J ' \n 2 
and to write the full noise covariance matrix as 
/ u ( Ni N12 \ 



(2) 



(3) 



We review the mathematical details of combining, fil- 
tering and comparing maps in Appendices A, B and C, 



with some explicit details added (beyond |28|]) that are 
useful when implementing these methods in practice. 

We derive a new deconvolution method in Appendix |d| 
which substantially simplifies our calculations by elimi- 
nating the A-matrices above. In the generic case, de- 
convolution is strictly speaking impossible, since the ma- 
trix A is not invertible and certain pieces of informa- 
tion about x are simply not present in y. It is common 
practice to find approximate solutions to such under- 
determined problems using singular value decomposition 
or other techniques, but our goal is different: we need a 
deconvolved sky map x that can be analyzed as a true 
sky map with A = I without approximations, shifting 
all complications into the new noise covariance matrix. 
The method derived in Appendix ^ is found to be sta- 
ble numerically, and can be used both for "unsmoothing" 
low-resolution maps and to deconvolve more complicated 
oscillatory beam patterns such as that of the Saskatoon 
experiment. 



III. RESULTS 

In this section, we combine the QMAP, Saskatoon and 
COBE maps. We then perform a battery of tests for sys- 
tematic errors by comparing the maps with each other, 
paying particular attention to possible calibration, point- 
ing and foreground problems. 



A. Saskatoon Data 

The Saskatoon data set is very different from other 
data sets such as QMAP and COBE since it does not con- 
tain simple sky temperature measurements. Instead, the 
2970 Saskatoon measurements are different linear combi- 
nations of the sky temperatures with rather complicated 
weight functions reminiscent of caterpillars — examples 
are plotted in |2l],|23|]. These measurements probe a cir- 
cular sky patch with about 16° diameter, centered on 
the the North Celestial Pole (NCP). In addition to the 
2590 weight functions used in Q , which are all oriented 
like spokes of a wheel, we include the 380 "RING" data 
measurements, linear combinations in the perpendicular 
direction going around the periphery of the observing re- 
gion. 

We pixelize this sky region into 2016 pixels in the same 
coordinate system as QMAP, i.e., a simple square grid 
in gnomonic equal area projection, and define x to be 
the true sky convolved with a Gaussian beam of FWHM 
0.68°. We compute the 2970 x 2016 matrix A using 
the software from the original Saskatoon analysis J3l],0 ■ 
The rows of A have a vanishing sum since the beam func- 
tions are all insensitive to the monopole — they are nor- 
malized so that the absolute values sum to two. Using 



equations (Dl) and (D2), we obtain a deconvolved Saska- 
toon CMB map xgASK and its corresponding covariance 
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matrix Xsask- Figure || shows the Wiener-filtered Saska- 
toon map, which is visually almost identical to that in 
p3| except for the additional information from the RING 
data near the border. 




FIG. 2. Wiener-filtered Saskatoon map. The CMB temperature 
is shown in coordinates where the north celestial pole is at the 
center of a circle of 16° diameter, with R.A. being zero at the 
top and increasing clockwise. In addition to the data included in 
the map of j^jj], "RING" data is included here. Note that the 
orientation of this and all following maps is different from that in 
Figure hi 




FIG. 3. Wiener-filtered QMAP map. The coordinates are the 
same as in the previous figure. 




FIG. 4. Wiener-filtered map combining the QMAP and Saska- 
toon experiments. The coordinates are the same as in Figure bl 




FIG. 5. Wiener-filtered map of the combined QMASK and 
COBE data. The coordinates are the same as in the previous fig- 
ure. COBE adds only large-scale information. For example, the 
upper left region is brightened somewhat. 
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B. Combining QMAP with Saskatoon 

The QMAP data consists of Ka-band and Q-band mea- 
surements with angular resolution 0.89° and 0.68°, re- 
spectively. We first deconvolve the Ka-band data to 0.68° 
resolution using the method of Appendix [D[ We then 
produce a unified QMAP data set (xqmapjSqmap) 
with a single resolution, 0.68°, by combining the result 
with the Q -band as described in Appendix |X|, using equa- 
tion (All) to compute the final noise covariance Hqmap 
since the overlap between the two bands is only partial. 
This combined QMAP map has 5396 pixels covering a 
sky area of about 538 square degrees. 

Combining QMAP and SASK is now a straightforward 
task, since the data sets xqmap and Xsask have the same 
angular resolution and pixelization scheme. Since the 
spatial overlap is only partial, we once again use equa- 
tion (All) to compute the combined noise covariance ma- 
trix. There are 917 overlapping pixels, so the combined 
map consists of 6495 pixels, covering a sky area of about 
648 square degrees. We will nickname the combined data 
set "QMASK" . 

The main improvement in the combined is not the 
area covered (the QMASK map is only 20% larger than 
QMAP), but the signal-to-noise and the topology. SASK 
has excellent signal-to-noise in the region that it covers, 
which overlaps the most sensitive region of QMAP. In- 
deed, the two maps have comparable sensitivity in the 
overlap region, so both of them have substantial impact 
on the spatial features seen in Fig. 3. Filling in the 
"hole" in the map is also useful for comparing with lower 
resolution maps like COBE and for potential future ap- 
plications, e.g., genus statistics, where a large contiguous 
area is desirable. 



C. Combining QMASK with COBE 

The COBE data (3^]34|] has much lower angular res- 
olution than QMASK (about 7.08°), and the pixel size 
of COBE is much bigger than that of QMASK as well 
(about 2.6° x 2.6°). In total, there are 6144 pixels in 
the whole COBE sky map. We select those pixels whose 
centers are within the QMASK map and at least 3° away 
from the perimeter. Only 58 COBE pixels satisfy these 
criteria. 

We first deconvolve this COBE data to the QMASK 
angular resolution, using the method of Appendix ^ with 
Ai being a 58 x 6495 matrix with a Gaussian COBE beam 
on each row, normalized to sum to unity. As input, we 
use the inverse-variance weighted average of the 53 and 
90 GHz COBE/DMR channels 0. By construction, our 
resulting COBE and QMASK maps overlap each other 
perfectly, s o we obtai n ou r combined map by simply using 
equations (A5) and (A6). 

At first glance, Figure [| and Figure || look very similar. 
However, inspecting them more carefully reveals that al- 



though the small scale patterns are the same, the upper 
part in Figure [| is brighter than that in Figure || (this 
is related t o the QMASK cold spot that we will discuss 
subsection [II E|) . In other words, since COBE contains 
only large scale information, this is precisely what it has 
added in Figure ^, leaving the small scale structure un- 
affected. 



D. Comparing QMAP with Saskatoon 

As mentioned above, there are 917 pixels overlapping 
between the QMAP and Saskatoon data sets. After ex- 
tracting these pixels and their two noise covariance ma- 
trices from the full maps, we can compare the two experi- 
ments using the null-buster test described in Appendix |^. 
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FIG. 6. Comparison of QMAP (left) and Saskatoon (right). 
The upper panel shows both maps Wiener filtered with the same 
weighting in the overlap region. The lower panel shows the num- 
ber of standard deviations ( "sigmas" ) at which the difference map 
xqmap — '"^sask is inconsistent with mere noise. Note that this 
is only for the overlap region. 
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1. Visual comparison 

Before delving into statistical details, it is useful to 
compare the two maps visually. Comparing plots of the 
raw maps x is rather useless, since they are so noisy. Un- 
fortunately, comparing plots of the two Wiener-filtered 
maps like figures ^[^| is not ideal either: since the noise 
matrices Ni and N2 are different, this would entail com- 
paring apples and oranges, since the two Wiener-filtered 
maps would be smoothed and weighted differently For 
instance, if a prominent spot in one map is invisible in 
the other, this could either signal a problem or be due 
to that particular region being very noisy in the second 
map and therefore suppressed by the Wiener filtering. 

To circumvent this problem, we Wiener-filter both 
maps exactly in the same way, 



x™ =S[S + N!+N 2 ] 



= 1,2. 



(4) 



In other words, we replace the individual noise covariance 
matrices by their sum, so that the map will only show 
information that is accurately measured by both experi- 
ments. These maps x™ and X2 are compared in the up- 
per part of Figure ^, and are seen to look encouragingly 
similar. 
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FIG. 7. Same as previous figure, but comparing Saskatoon with 
various subsets of the QMAP data. The curves correspond to 
Ka-band (dotted), Q-band (solid), flight 1 (short-dashed), flight 2 
(long-dashed), f2Kal2(dot-long-dashed), f2Q12 (dot-short-dashed) 
and f2Q34 (short-dashcd-long dashed). 



2. Tests for systematic and calibration errors 

The lower plot in Figure || shows the results of apply- 
ing the null-buster test to the difference map xqmap — 
rxgASK for different values of the constant r. (The corre- 
sponding noise covariance matrix N = Ni + r 2 N2.) The 
left part of the curve where r <C 1 is dominated by infor- 
mation from QMAP, and we see that QMAP alone (the 



r = case) is inconsistent with noise at about the 40u- 
level. Similarly, we see that the Saskatoon map alone (the 
case r = 00) is inconsistent with noise at about the 20a- 
level. Note that these significance levels are still higher 
for the full maps — here we are limiting ourselves to the 
sky region where they overlap. 

In summary, both the QMAP and Saskatoon maps 
contain plenty of signal. Is this signal consistent between 
the two maps? The answer is given by the most interest- 
ing point on the curve, where the relative normalization 
value r = 1. In the absence of systematic or calibration 
errors, the corresponding difference map should contain 
pure noise. In our case, when r = 1, the difference map 
xqmap — xsask is seen to be consistent with pure noise, 
i.e., less than 2a away from zero. The strong signal seen 
in both maps therefore appears to be a true sky signal, 
with no evidence for significant systematic errors in ei- 
ther QMAP or Saskatoon. 

The QMAP experiment consists of two flights, each 
with two frequencies (Ka and Q-band) and three slightly 
different observing regions. We label these six sub-maps 
flKal2, flQ2 and flQ34 (from flight 1) and f2Kal2, 
f2Q12 and f2Q34 (from flight 2). The pointing and cali- 
bration analyses for these two flights were completely sep- 
arate, and the map-making algorithm was applied sepa- 
rately for these six sub-maps |25|-|27j . To investigate pos- 
sible problems with these individual sub-maps that may 
have been averaged away in the combined analysis, we re- 
peat the comparison with Saskatoon separately for each 
one. The results are shown in Figure [?]■ None of these 
curves show any evidence for systematic or calibration 
errors. The f2Q34 map is seen to contain the strongest 
signal, inconsistent with noise at the 38cr level. This is 
because f2Q34 contains a striking cold spot — we will 



return to this in more detail in subsection III G 



Table 1. Summary of map comparisons. The first three lines 
give the number of "sigmas" at which map 1, the difference map 
and map 2, respectively, are inconsistent with noise. The remaining 
lines give the best fit value and limits on the relative calibration r, 
or, for the Ka vs. Q case, the spectral index j3. 





QMAP vs SK 


QMASK vs COBE 


Ka vs Q band 


v(r = 0) 


40 


62 


40 


v\t = 1) 


1.97 


-0.64 


-0.53 


u(r = 00) 


21 


3 


26 


T min 


0.79 


0.09 


-2.73* 


T min (2c) 


0.95 


0.13 


-2.2* 


Train (if ) 




0.2 


-1.7* 


Tbest 


1.1 


1.25 


0.* 


T max (if ) 




11.7 


1.42* 


Tmax (2<t) 


1.2 


20 


1.9* 


Tmax (3<t) 


1.48 


63 


2.2* 



* values of the power spectral index j3 in subsection 111 F 



The null-buster curves are interesting at more than just 
the points r = 0, 1, and 00; the entire region near r = 1 
places limits on calibration errors. A relative calibration 
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error of say 10% would shift the minimum of the curve 
sideways to 0.9 or 1.1, depending on whether the QMAP 
or Saskatoon map was too high. We can therefore place 
limits on calibration errors by reading off the r-values 
where noise is ruled out at say 2 sigma. For instance, 
the QMAP-Saskatoon comparison constrains the relative 
calibration error to be less than 20% (3c), and the re- 
values for different significance levels are shown in table 
1. This method may prove quite useful for upcoming 
experiments that have higher sensitivity. 

3. Pointing tests 

Our comparison method can also be used to test for 
relative pointing errors, as a complement to the standard 
lower-level pointing tests that are routinely made using 
point sources etc. Although an overall sideways shift of a 
single map will not affect the measured power spectrum, 
such errors can become disastrous if the map is combined 
with another one. 

As an illustration of such a test, we compare the f2Q34 
map with the Saskatoon map with the null-buster test 
(setting r = 1) after shifting it vertically and horizontally 
by an integer number of pixels. Figure @ shows that there 
is no evidence for pointing error although we cannot give 
a strong constraint. Just as the calibration test, this 
pointing test based on CMB maps alone is likely to be 
useful for upcoming high-sensitivity experiments. 
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FIG. 8. Test of the relative pointing of QMAP and Saskatoon. 
The curves show the number of "sigmas" at which the difference 
map is inconsistent with noise when the QMAP map is shifted 
vertically and horizontally. Starting from the inside, the contours 
are at 1, 2, 3, 4, 5, and 6cr respectively. Cross indicates no shift. 
The pixels are squares of side 0.3125°. 



FIG. 9. Comparison of QMASK (left) with COBE (right). The 
upper panel shows both maps Wiener filtered with the same weight- 
ing in the overlap region, using equation (Q). The lower panel shows 
the number of standard deviations ( "sigmas" ) at which the differ- 
ence map xqmask — rxcoBE is inconsistent with mere noise, and 
illustrates that the visual discrepancy at "2 o'clock" is consistent 
with a fluctuation in the (correlated) noise. 



E. Comparing QMASK with COBE 

Our QMASK map and our deconvolved COBE map 
cover the exact same sky region, so there are 6495 pix- 
els in the overlap maps. Generating Figure ^|, which 
compares these two maps, therefore involved a marathon 
computer run, processing the 6495 x 6495 matrices of 
equation (CI) for each r-value.^] As can be seen in the 
lower panel of Figure ||, the QMASK data is inconsistent 



* If CPU time had been an issue, this calculation could have 
been accelerated by binning the QMASK pixels into larger 
ones. This would give essentially the same answer, since the 
null-buster test gives statistical weight only to modes where 
both maps are sensitive — in this case, to large scale modes 
only. 
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with noise at about the 62c-level, whereas the COBE 
data is inconsistent with noise slightly above the 3cr- 
level. Since neither the QMAP nor the Saskatoon experi- 
ments were designed to probe such large angular scales, it 
is quite encouraging that the QMASK-COBE difference 
map (r = 1) is seen to be consistent with pure noise. 
Since the minimum of the curve is so broad, however, we 
obtain no interesting constraints on calibration errors. 

The upper panel in Figure |^ shows that the two maps 
look fairly similar considering the weak (3er) COBE sig- 
nal, with the notable exception of the upper right part 
of the Saskatoon disk. Here COBE shows a hot spot 
whereas QMASK shows a cold spot. We will re turn to 
this issue in more detail below, in subsection [II G . When 



generating these two maps, we used the same equal- 
weighting Wiene r filteri ng method that was described 
in subsubsection HID 1 . This is particularly important 
here, since the QMASK and COBE have such dramati- 
cally different angular resolutions — in contrast, a visual 
comparison of the Wiener filtered COBE map with the 
normal Wiener filtered QMASK map from Figure [| is 
rather useless, since the latter is dominated by small- 
scale fluctuations. 



F. Foreground constraints 

The previous two sections used map comparisons to 
test for calibration and pointing errors. Here we will 
compare maps at different frequencies to constrain the 
spectrum of the detected sky signal. 

The presence of foreground contamination (see |35| l 
for a recent review) has been quantified for both the 
Saskatoon [||| and QMAP [|6| experiments by cross- 
correlating the maps with various foreground templates. 
The dominant foreground emission is expected to be due 
to synchrotron radiation, free-free emission and (vibra- 
tional and spinning) dust emission, from both the Milky 
Way (seen as diffuse emission) and other galaxies (seen 
as point sources). These cross-correlation analyses con- 
cluded that foregrounds played only a subdominant role 
in Saskatoon and QMAP. 

By comparing the Ka- and Q-band maps, we are able 
to place a direct constraint on the frequency dependence 
of the signal. Both Saskatoon and QMAP observe in 
both of these frequency bands. We therefore repeat 
the analysis described above (Saskatoon deconvolution, 
merging with QMAP, etc.) separately for each of the 
two bands. The upper panel of Figure [lO] shows equal- 
weighting Wiener-filtered maps for Ka-band and Q-band 
in the sky region that was observed at both frequencies, 
showing that they visually agree well. 

If we fit the frequency dependence by a power law 
8T(v) (x v@ over the narrow frequency range in ques- 
tion (fKa ~ 30 GHz, i/Ka ~ 40 GHz), then characteristic 
spectral indices are f3 ~ —2.8 for synchrotron, (3 ~ —2.15 
for free- free emission, j3 ~ — 3 for spinning dust and (3^2 



for vibrating dust. By definition, (3 = for CMB. If the 
sky signal in our maps obeyed 5T{v) cx , then the dif- 
ference map XKa — ricQ would contain pure noise when r 
was such that 



Ig^Ka/^Q) ' 



(5) 



The lower panel of Figure [Tc] is therefore plotted with j3 
rather than r on the horizontal axis. Insisting that the 
difference map not be inconsistent with noise at more 
than lcr gives the spectral index constraint = O.Ol^'y, 
which is inconsistent with the signal being any single one 
of the foregrounds mentioned above. 
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FIG. 10. Comparison of QMASK at two different frequencies, 
Ka-band (left) and Q-band (right). The upper panel shows both 
maps Wiener filtered with the same weighting in the overlap region 
that was observed at both frequencies. Arrows indicate the coldest 
spot discussed below. The lower panel shows the number of stan- 
dard deviations ("sigmas") at which the difference map XK a - rxq 
is inconsistent with mere noise. 



G. The coldest spot 

Up until now, we have presented a battery of tests for 
systematic errors and other problems, all of which passed. 
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However, our maps did turn up one somewhat anomalous 
feature: an unusually cold spot around "two o'clock" 
in the Saskatoon disk. The spot's coldest pixel in the 
QMASK map is located at RA=3 /l 20 m , DEC=84°55'. 
Here the Wiener-filtered Q-band map in Figure |l0| gives 
ST w — 230/iK. For comparison, the expected rms fluc- 
tuations in this map are 27/iK from detector noise and 
49/iK from CMB fluctuations (for the "concordance" 
power spectrum of pJ|), summing to 56/xK in quadra- 
ture. Taken at face value, this would indicate that the 
spot is a —4.1a fluctuation. For completeness, this sec- 
tion describes a number of additional tests performed in 
an attempt to clarify its nature. 

It is unlikely that the cold spot is due entirely to sys- 
tematic error, since it is clearly detected by both QMAP 
(covered by the Q3 and Q4 detectors from in flight 2) and 
Saskatoon. The Saskatoon experiment even detected this 
spot independently in each of its three observing seasons 
pl| , p3[ . Indeed, the reason that the f2Q34 map shows the 
most spectacular agreement with Saskatoon in Figure 
is that f2Q34 is the only QMAP map that covers the area 
containing this spot. 
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FIG. 11. The temperature towards RA=3 h 20 m , DEC=84°55' 
at different frequencies. Errors bars correspond to detector noise 
alone. Note that these points cannot be interpreted as a spectrum 
of this sky region, since the COBE points have much lower reso- 
lution and the two QMASK points have been Wiener filtered with 
different weights, pushing them closer to zero than the underlying 
sky temperature. 

We have plotted all available microwave measurements 
of this region in Figure ^ as a function of frequency, in- 
cluding both the Ka- and Q-band measurements from 
QMASK and the COBE/DMR observations at 31.5, 53 
and 90 GHz |l7j]. Unfortunately, Figure |ll| is not an 
actual spectrum of the spot. It is more of a compari- 
son of apples and oranges, since the measurements differ 
dramatically in angular resolution (7° for COBE) and 
the QMASK maps are Wiener-filtered. Since Wiener fil- 
tering always pushes the signal towards zero when noise 
is present, the QMASK points should be interpreted as 



lower limits on |<5T| — the true sky temperature is likely 
to be even colder. 

On the seven degree scale probed by COBE, no evi- 
dence is seen for a cold spot (this can also be seen in the 
map of Figure ||) , and the COBE spectrum of this region 
appears consistent with CMB, i.e., (3 = 0. One pos- 
sible interpretation is that a small cold region resolved 
by QMAP and Saskatoon is partly smoothed away by 
COBE. 

Unfortunately, this sky patch is not covered by QMAP 
in Ka-band, so the Ka-information comes from Saska- 
toon alone and is therefore noisier than the Q-band mea- 
surement. This means that the Wiener-filtering has sup- 
pressed the Ka-band more, as well as lowered its res- 
olution by more aggressive smoothing. However, the 
expected rms CMB fluctuations in the Wiener-filtered 
maps are not nearly as different as the data points in 
Figure |l"l| (— 31^K at Ka-band and —230/iK at Q-band), 
indicating that the low Q-band temperature does not per- 
sist fully down to Ka-band. This argues against both a 
CMB origin and a thermal SZ-origin, which would cause 
a cold spot that was essentially frequency-independent 
for v < 100 GHz. 

All other known microwave foregrounds produce hot 
rather than cold spots. Some sort of absorption process 
also appears unlikely, since the absorbing medium would 
have to be colder than 3K. No relevant foreground emis- 
sion or X-ray cluster is found in radio, infrared or X-ray 
maps of the region. Some adjacent dust emission is seen 
in the IRAS 100/im map |32[| , which could potentially 
make this region look cold in contrast (since none of our 
maps are sensitive to the monopole mode and measure 
merely relative temperatures), but only at the level of 
about 10 /iK. 

In Q-band, the spot is observed by QMAP only in 
flight two and only in the Q3 and Q4 channels. The lat- 
ter dominates statistically, and has a 13% calibration un- 
certainty. However, Saskatoon also observes the spot in 
Q-band, and both experiments measure around — 200/xK 
in their individual Wiener filtered maps. 

The explanation is fairly likely to be something mun- 
dane, since a 4-sigma fluctuation (which should hap- 
pen about once for every 16,000 independent pixels) is 
not extremely unlikely when many pixels are consid- 
ered — if the effective number of independent regions 
in the Wiener filtered map is 10 2 taking the smoothing 
into account, the significance level gets downgraded to 
(1 - 1/16000) 100 S3 99% |3 • 

In summary, the cold spot is definitely out there at 
some level, but we have no single simple interpretation 
of what is causing it. Its unusual spectrum argues against 
a non-Gaussian CMB fluctuation, foreground contamina- 
tion and an SZ signal. The most likely remaining expla- 
nation is a confluence of a less extreme CMB cold spot, 
noise fluctuations, calibration uncertainty and perhaps 
some small systematic error. We have described this spot 
in such detail simply to ensure that no hints of problems 
with the data get swept under the rug. 
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The MAP satellite should resolve this puzzle next year, 
observing the spot with high sensitivity and resolution at 
22, 30, 40, 60 and 90 GHz. 



IV. DISCUSSION 

We have presented methods for comparing and com- 
bining CMB maps and applied them to the QMAP, 
Saskatoon and COBE DMR data sets. We found that 
these methods were able to place interesting constraints 
on calibration and pointing problems, foreground con- 
taminations and systematic errors in general. This 
should make them quite useful for ongoing and upcoming 
high-precision experiments. 

The data sets passed our entire battery of consistency 
tests, placing strong limits on systematic errors. Al- 
though the Saskatoon and QMAP maps detect signal at 
the 21<7 and 40ct levels in the overlap region, respectively, 
their difference is consistent with pure noise. Our com- 
bined QMAP + Saskatoon map, nicknamed QMASK, 
covers a large enough area to allow a statistical compar- 
ison with COBE/DMR, showing good agreement. 

The one surprise that our battery of tests turned up 
is that a small region around (RA=3 h 20 m ,DEC=84°55') 
appears unusually cold, mainly in Q-band. Its unusual 
frequency dependence argues against non-Gaussian CMB 
fluctuations, SZ-signal and known foregrounds. 

The QMASK map presented here has been made pub- 
licly available at www. hep. upenn.edu/^xuyz/qmask. html 
together with its 6495 x 6495 noise covariance matrix. 
With its 648 square degrees, this thoroughly tested data 
set is currently the largest degree-scale map available, 
detecting signal at a level exceeding 60a. 
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6034; Hubble Fellowships HF-01044.01-93A (to TH) and 
HF-01084.01-96A (to MT); a Cottrell Award from the 
Research Corporation (to LP) and a David and Lucile 
Packard Foundation Fellowship (to LP). The COBE data 
sets were developed by the NASA Goddard Space Flight 
Center under the guidance of the COBE Science Working 
Group and were provided by the NSSDC. 



APPENDIX A: COMBINING MAPS 



The combined map x defined by 

x= [A*N _1 A] _1 A'N _1 y 



(Al) 



can be shown to be unbiased ((x) = x), to minimize the 
rms noise in each pixel and, if the noise properties are 
Gaussian, to retain all information about the true sky 
x that was present in the two original maps |28f| . The 
corresponding covariance matrix of the noise e = x — x is 



(ee 4 ) = [A'N^A]" 



(A2) 



In all cases treated in this paper, the noise is uncorrelated 
between different maps ((N12 = 0), which simplifies these 
equations to 



x = S [AfNfVi + AIN2V2] 
£ = [ASNJ^Ai + AlN^Aa] " 



(A3) 
(A4) 



Thanks to the deconvolution technique that will be de- 
scribed in Appendix we will generally face the much 
simpler case where the two data sets are two sky maps 
with the exact same angular resolutions, i.e., the case 
Ai = A2 = I, reducing the last two equations to simply 



x = s [NrVi + N2V2] 

S= [N^+Nj 1 ] -1 . 



(A5) 
(A6) 



For the case of only a single pixel, we recognize this as a 
familiar inverse- variance weighting. More generally, if the 
two noise matrices can be simultaneously diagonalized, 
we see that this combination scheme corresponds to an 
inverse-variance weighting eigenmode by eigenmodc. 

Generally maps overlap only partially, so we need only 
apply this matrix method in the common region. Yet 
care needs to be taken in computing the noise covariance 
matrix S of the final map, since it will contain correla- 
tions between the common region and the rest. Let us 
split the noise vectors for the two maps as 



ni = 



n c i 



n 2 



n c2 



(A7) 



where the subscript c refers to the common region of 
the two maps whereas a and b refer to the regions that 
only belong to maps 1 and 2, respectively. We write the 
corresponding covariance matrices as 



Ni 



No 



{n 2 n 2 ) 



( n el n a) 



N cl 



n c2 



(A8) 



(A9) 



Substituting equation ([y) into equation ([A3]), we obtain 
the noise vector n c for the combined map in the common 
region: 



n c = S c [A^N^n^ + A* c2 N^n c2 ] 



(A10) 



where S c is given by equation ( A4) for the common part. 
The final combined noise covariance matrix S for the 
noise vector (n Q , nt,, n c ) of the combined map is therefore 
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/ N„ (nan*>\ 

S = N b (n b n*) , (All) 

UW) ([n 6 n*] 4 ) S c 7 

where 

(n a n* c ) = (n Q n* 1 )N- 1 1 A cl S c , (A12) 
(n b n^> = (n^N^A^E,. (A13) 

We need to use these expressions repeatedly in this paper 
to combine partially overlapping maps, e.g., combining 
the different QMAP flights with each other and combin- 
ing QMAP with Saskatoon. 

APPENDIX B: PLOTTING MAPS 

Although the map x contains all the sky information 
from an experiment, plotting it is not very useful when 
some modes are much more noisy than others, thereby 
dominating the visual image. For this reason, it has be- 
come standard in the community to plot the correspond- 
ing Wiener filtered map, defined as 

x lu = SIS + N]- 1 *, (Bl) 

where S is an estimate of the covariance matrix due to sky 
signal. Throughout this paper, we use the S-matrix cor- 
responding to the "concordance" power spectrum from 
PH ) which agrees well with all current CMB measure- 
ments. 



be interpreted as the number of "sigmas" at which Hq is 
ruled out [p8| . Note that for the special case S a N, it 
reduces to simply v = (x 2 — n)/\/2n, where x 2 = z*N _1 z 
is a standard chi-squared statistic. The null-buster test 
can therefore be viewed as a generalized x 2 ~test which 
places more weight on those particular modes where the 
expected signal-to-noise is high. It has proven successful 
comparing both microwave background maps p5|-p7j] and 
galaxy distribution [ p9|j30|. 

To evaluate equation( |Cl| ) in practice, it is useful to 
Cholesky decompose the noise matrix as N = LL* and 
compute the matrix R = L _1 SL _ *. The remainder of 
the calculation now becomes trivial, since tr{N _1 S} = 
trR = £R 4l and trjN^SN^S} = trR 2 = £(R 4j ) 2 . 

For the general case when Ai ^ A2, the situation 
is more complicated, since it is non-trivial to construct 
a difference map which is free of sky signal. A tech- 
nique involving signal-to-noise eigenmode analysis has 
been derived for this case (2Sj| , but it is unfortunately 
rather complicated and cumbersome to implement. Be- 
low we present a simpler method that eliminates need for 
this by reducing the general problem to the simple case 
A 1= A 2 =I. 

APPENDIX D: DECONVOLVING MAPS 

In this section, we present a method for inverting equa- 
tion ([!]), i.e., for undoing the convolution with beam and 
scanning effects given by the A-matrices. 



APPENDIX C: COMPARING MAPS 



Why is it useful? 



Here we discuss the issue of how to test whether two 
data sets are consistent or display evidence of systematic 
errors. Specifically, is there some true sky x such that 
the data sets yi and y 2 are consistent with equation |l])? 
Let us first consider the simplest case where the two data 
sets sample the sky in the same way, that is, Ai = A 2 . 
Consider two hypotheses: 

Ho: The null hypothesis Ho that there are no system- 
atic errors, so that the difference map z = yi — y 2 
consists of pure noise with zero mean and covari- 
ance matrix (zz*) = N = Ni + N 2 . 

H±: The alternative hypothesis that the difference map 
z consists of some signal besides noise, i.e., (z) = 0, 
and (zz 4 ) = N+S for some signal covariance matrix 
S. 

The "null-buster" statistic |E1 



z'N-iSN^z - tr {N^S} 



[2 tr{N- 1 SN- 1 S}] 



1/2 



(CI) 



can be shown to rule out the null hypothesis Ho with 
the largest average significance {v) if Hi is true, and can 



As we will see, this simplifies calculations by eliminat- 
ing all A-matrices, encoding the corresponding compli- 
cations and correlations in the noise covariance matrices. 

When comparing or combining two maps, it is gen- 
erally undesirable to smooth the higher resolution one 
down to the lower resolution of the other, since this de- 
stroys information. Moreover, this tends to cause numer- 
ical instabilities by making the smoothed noise covariance 
matrix poorly conditioned. We will see that, surprisingly, 
deconvolution can be better conditioned than convolu- 
tion/smoothing. 

This deconvolution (elimination of A-matrices) is use- 
ful not only for comparing data sets as mentioned above, 
but for combining them as well. The complication stems 
from the fact that the sky is sampled by only a finite 
number of pixels, so to avoid problems with undersam- 
pling, x must be the "true sky" beam-smoothed map with 
some finite angular resolution. If we do not deconvolve, 



but use equations (A3) and (A4) to combine two data 
sets with different angular resolutions, say B\ = 0.89° 
and 6 2 = 0.68° for the QMAP Ka- and Q-band maps, 
respectively, it is natural to define x to have the higher 
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(0.68°) resolution, setting A 2 = I and letting Ai incor- 
porate the extra smoothing A6> = y/0\ — 9\ in the lower 
resolution map. The resulting map will now contain two 
kinds of pixels: ones with resolution # 2 in the region cov- 
ered only by map 2 and with the higher resolution 9\ 
elsewhere. If we need to combine this with a third map, 
the relevant smoothing scale unfortunately becomes un- 
defined near the boundary between the two resolutions. 
Deconvolution eliminates all these problems. 



b. How does it work? 

In the generic case, deconvolution is strictly speaking 
impossible, since the matrix A is not invertiblejj] and cer- 
tain pieces of information about x are simply not present 
in y. It is common practice to find approximate solutions 
to such under-determined problems using singular value 
decomposition or other techniques, but our goal here is 
different. We wish to compute a vector x that can be 
analyzed as a true sky map. Specifically, we want analy- 
sis of (x, S) to give exactly the same results as analysis 
of (y, N, A) for all cosmological applications, say Wiener 
filtering or power spectrum estimation. 

Our basic idea is to accept that certain modes in the 
map x cannot be recovered, and to record this infor- 
mation in the noise covariance matrix S by assigning 
a huge variance to these modes. Any subsequent analy- 
sis (say Wiener filtering or power spectrum estimation) 
will then automatically assign essentially zero weight to 
these modes. 

In practice, we find it conceptually useful to imagine 
combining our data yi with a "virtual map" y 2 that is 
so noisy that it contains essentially no information, yet 
has the angular resolution 9 that we wish to deconvolve 
down to, i.e., A 2 = I. Although, as we will see, this 
virtual map never enters the calculations in practice, it 
is a useful notion for intuitively understanding what the 
deconvolution technique does. Specifically, let us take the 
noise in the virtual map to be uncorrelated, with noise 
covariance matrix N 2 = u 2 l for some very large noise 
level a. Equations ( |A3| ) and (A4) give 



x = SA' 1 Nr 1 yi+^ 2 Sy 2 , 
X = [AiN^Ai + a~ 2 l] -1 . 



(Dl) 
(D2) 



In the limit a > 00, x will clearly become independent 
of the virtual temperature map y 2 except for the "junk 
modes" which have infinite variance according to X. For 
convenience, we therefore set y 2 = in practice.^] 



' Specifically, the problem is that A generically has a non- 
zero null space, i.e., that there are non-zero vectors x such 
that Ax = 0. 

'''An alternative approach would be to set V2 equal to a 



This deconvolution method has exactly the property 
we want as long as a is orders of magnitude larger than 
the pixel signal due to CMB. If we were to choose a to be 
too small, then the virtual map would contribute a non- 
negligible amount of information and bias the results. If 
we were to choose a to be too large, however, the ma- 
trix X would contain some enormous eigenvalues (since 
A^N]Ai is typically not invertible) and be poorly condi- 
tioned, which could cause numerical problems in subse- 
quent analysis. We performed a series of numerical tests 
to assess these problems, and found that with n < 10 4 
pixels and double precision arithmetic, a = 10 4 /iK was 
a good compromise that produced neither of these two 
problems. We will therefore use this choice throughout 
the present paper. 



c. Tests 

As a first test of the method, we deconvolve (or "un- 
smooth") a map yi with resolution 9\ into a map with 
resolution 9 2 (82 < 61). For this case, 



(Ai), 



1 



2ttA6> 2 



(D3) 



where 



cos 1 (?i • Vj ) is the angular separation be- 



tween pixels i and j, and A9 = \/9\ — 9\ is the extra 
smoothing to be undone. Specifically, we unsmooth the 
QMAP Ka-band data to obtain the same resolution as 
the Q-band data has, from 0.89° to 0.68°. We then 
Wiener filtered both the original and unsmoothed ver- 
sions of the map, obtaining virtually identical results. 

As a second test, we deconvolve the Saskatoon data 
y into a map x using the full (and rather complicated) 
A-matrix described in ||J . We then Wiener-filter x and 
obtain a map virtually identical to the one that was com- 
puted in J23j — the latter was computed with a com- 
pletely different method which circumvented the map 
step altogether. The details of these maps have already 
been presented above. 

Both of these tests thus confirm what we expect the- 
oretically: that the deconvolved map (x, S) contains ex- 
actly the same information about the true sky as the 
input data (y,N, A), no more and no less. 



Monte-Carlo generated map of independent Gaussian random 
variables with standard deviation a — although this results 
in different numerical values in x, it will of course not change 
the results of any subsequent cosmological analysis of x, since 
only the "junk modes" are different. 
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